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Abstract 

A new post-processing technique has been devel- 
oped, based on the Intrinsic Low Dimensional Man- 
ifold (ILDM) method of Maas and Pope 1,2 . 

The ILDM method is a dynamical systems ap- 
proach to the simplificat ion of large chemical kinetic 
mechanisms. By identifying low-dimensional at- 
tracting manifolds, the method allows complex full 
mechanisms to be parameterized by just, a few vari- 
ables: In effect, generating reduced chemical mech- 
anisms by an automatic procedure. These resulting 
mechanisms however, still retain all the species used 
in the full mechanism. 

The NOx post-processor takes an ILDM reduced 
mechanism and attempts to map this mechanism 
to the results of a CFD calculation. This mapping 
allows the NOx concentrations at each grid node 
to be obtained from the ILDM reduced mechanism, 
as well as other trace species of interest. Because a 
mapping procedure is used, this method is very fast, 
being able to process one million node calculations 
in just a few minutes. 

Introduction 

The ability to predict NOx in CFD calculations 
of reacting flows is constrained by the trade-off be- 
tween accuracy and available computing resources. 
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For fluid design calculations, it is enough to get 
the values of density and temperature correct, and 
so fairly simple chemical mechanisms may be em- 
ployed, with good results. 

However the prediction of trace species, such as 
NO, requires the use of a significantly more complex 
kinetics scheme, with the resulting increase in com- 
putational effort required to solve the problem. For 
example, a simple two-step mechanism will have six 
species (Fuel, H20, CO, C02, 02, N2) and two rate 
equations. To predict NO by the Zeldovich mecha- 
nism will require adding 5 more species (OH O H N 
NO) and six more rate equations. Clearly this has 
more than doubled the complexity of the problem 
and the CPU time required. 

To overcome this increase in complexity, post- 
processing is an attractive option, with the promise 
of quickly obtaining the NOx fields when the so- 
lution of the major species and velocity fields has 
converged. To do this, two main approaches are 
employed. 

In the first, the composition and temperature at 
each node are used to obtain the equilibrium value 
of NOx. This has the advantage of being very 
quick, but will naturally over-predict the value of 
NOx at each point. Due to the very slow rate of 
NOx formation under most conditions, it may sig- 
nificantly over-predict the value. The advantage of 
this method is its speed, and that it does not need 
to know anything about the composition at neig- 
bouring nodes. 

The second technique employed is an iterative, 
frozen-composition method. In this the values of 
velocity, density and major species compositions 
are frozen and the species transport equations are 


NASA CR-204137 


1 



solved for the minor species. The disadvantage of 
this technique is the time required to perform the 
iterat ions to convergence. Typically this can be of a 
similar CPU time as the original calculation. In ad- 
dition, the technique is not code independent, i.e. 
The CFD results from an unstructured code will 
have to be processed differently than those from a 
structured code. The advantage is that the NOx 
reaction rates are not assumed to be fast, and so 
finite rate effects are accounted for. 

In this paper, an approach is outlined that at- 
tempts to provide the equilibrium approach with 
the finite rate effects available to the iterative 
scheme. 

To achieve this, an Intrinsic Low Dimensional 
Manifold (ILDM) simplificat ion of a full mechanism 
is performed, so that all species are represented as 
functions of a mixture fraction £/ and a progress 
variable Ypj . Then for each node in the CFD data, 
a mixture fraction and progress variable Ypc is 
obtained and by mapping the ILDM variable onto 
the CFD values, a value for the trace species may 
be obtained. 


ILDM Method 

The ILDM method employed in this work is that 
used in the NASA Lewis ILDM code. 3 In this code, 
the ILDM simplified mechanism is obtained by a 
trajectory-generated technique, described by Pope 
and Maas 4 . In this method, the full mechanism 
is parameterized by two scalars, a mixture fraction 
based on atom concentrations and a progress vari- 
able, based on either species mass fractions, tem- 
perature or Gibbs function. The resulting species 
concentrations, rates and properties are stored in 
look-up tables. A more detailed description can be 
obtained by Norris 3 . 

One concern in the creation of the ILDM tables 
is that because the reaction rate of NOx is very 
slow, the ability to interpolate accurately may be 
compromised near the fully burnt limit of Ypj. This 
is because the progress variable will only change a 
small amount, while the NOx values can change by 
an order of magnit ude. To counter this, the look-up 
table has been modified to cluster the nodes near 
the fully burnt value. In addition, the choice of 
temperature as the progress variable is found to give 
the best resolution for NOx values. 


Post-Processor 

The post-processor performs the functions of 
reading in the CFD data, obtaining the values of 
£c and Ypc , and then mapping the composition at 
each node to the manifold. 

The mixture fraction of the CFD data can be 
obtained by several methods, each based on atom 
population. The mixture fraction based on the m 
atom, £ m is given by 

Y - Y 

£ 1 m ± m } air ( t ^ 

— T7 —r r - ( 1 ) 

*m,fuel * m,air 

where Y m is the mass fraction of the m atom, 
Ymjuel is the mass fraction of the m atom in the 
fuel stream and Y m j U€ i is the mass fraction of the 
m atom in the air stream. 

The three options for progress variable are 
species, temperature or Gibbs function based. The 
species option allows one to obtain the progress 
variable as a normalized mass fraction of a chemical 
species, or the sum of several chemical species. For 
example one could chose a mixture fraction based 
on the mass fraction of C02, or the sum of the mass 
fractions of C02 and H20. The only restriction on 
this choice is that the resulting quantity should be 
single-valued. For this reason, a progress variable 
based on the mass fraction of OH would not be 
acceptable, as the OH value rises, peaks and then 
drops down to an equilibrium value. Another op- 
tion is to use temperature as a progress variable, 
while the Gibbs function provides the third option. 
In all of these three cases, the progress variable 
must be the same as the one used to obtain the 
ILDM table. 

Having obtained £c and Ypc for each CFD data 
point, the NOx value is simply obtained by inter- 
polating from the ILDM table at the equivalent £ j 
and Ypj locations. 

Test Cases 

Three test cases were employed to evaluate the 
performance of the post-processor: A hydrogen- 
air can combustor, a methane-air combustor and 
a jeta-air problem. 

The hydrogen-air combustor consists of a stepped 
tube, with swirling air being injected in radially. 
Hydrogen fuel is injected in downstream of the air 
inlet, without swirl. Initial conditions for this case 
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are air at 804K, hydrogen at 305K and a pressure 
of 20 atm. In addition, several other inlet condi- 
tions were investigated, corrosponding to different 
operation conditions. 

Due to the sensitive nature of the methane and 
Jet A cases, the configuration of these problems will 
not be reported. However the initial conditions for 
the methane case are a fuel inlet temperature of 
305K, an air temperature of 950K and a pressure of 
20 atm. The Jet A case has initial air temperat ures 
of 740K, with a fuel temperature of 350K and a 
pressure of 10 atm. 

Data available for these flows consists of exper- 
imental measurements of NOx levels at the exit 
plane of the hydrogen-air cases, while the methane 
and Jet A configurations have the values of NOx ob- 
tained at each data point by the inclusion of the 
Zeldovich mechanism in their reduced mechanism. 
It should be noted that the CFD results for all but 
the first hydrogen-air setup were not run to a fully 
converged condition. However the results are in- 
cluded for completeness. 

The full mechanisms used to obtain the ILDM ta- 
bles are the hydrogen-air mechanism of Miller and 
Bowman 5 , the methane-air mechanism of Miller 
and Bowman 5 and the JetA-air mechanism of Koll- 
rack 6 . 

Results 

To address the accuracy of the mapping proce- 
dure, the scalar fields for known major species were 
extracted by the postprocessor and compared to the 
CFD results. The mixture fraction was obtained 
based on the mass fraction of the H atom for the 
H2-air case, and the C atom for the methane and 
jetA cases. In all cases, the temperature was used 
as the progress variable. 

For the hydrogen-air Case 1, the predicted mass 
fraction of H20 is plotted against the CFD results 
in Fig 1. As can be seen, there is a good agreement 
between the CFD result and the manifold predic- 
tion. For the methane case, the predicted mass frac- 
tions of C02 and H20 are compared to the CFD 
results in Fig 2. and 3. respectively. Again there is 
very good agreement, except for the very low val- 
ues. It can be seen that the values of C02 are 
underpredicted while the mass fractions of H20 are 
overpredicted. For the JetA example, the predicted 
mass fractions of C02 and H20 are compared to 
the CFD results in Fig 4. and 5. respectively. In 
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this case, the values of C02 are underpredicted for 
the entire range of values, while the H20 values are 
overpredicted. This deviation in the results will be 
discussed in the next section. 

For the hydrogen-air cases, the NO exit concen- 
tration was predicted and compared to experimen- 
tal data. The results of this are shown in Table 1. It 
can be seen that a reasonable agreement is obtained 
in all 5 cases. It should be noted that the first case 
is the only one to have converged fully, thus the re- 
sults for the other four may be somewhat suspect. 
For the hydrogen case, no CFD values of NOx were 
available. 

Hydrogen-Air NO Prediction Performance. 


Case Number 

Prediction 

Experiment 

Case 1 

20.5 ppm 

21.6-26.5 ppm 

Case 2 

6.9 ppm 

10.6 ppm 

Case 3 

4.9 ppm 

12.0 ppm 

Case 4 

13.9 ppm 

3.1 ppm 

Case 5 

28.5 ppm 

35.1 ppm 


For the methane and JetA case, there are no ex- 
perimental results available, and so the NOx pre- 
dicted by the post-processor was compared to the 
CFD values. 

In the Methane case, a reasonable agreement was 
obtained between the the prediction and the CFD 
results. A small percentage of results were under- 
predicted, and can be seen as the band of points 
lying to the left of the diagonal line in Fig 6. 

For the JetA case shown in Fig. 7, the predicted 
NOx does not agree well with the CFD data, except 
at the high concentration levels. 

Discussion 

The results described in the previous section will 
now be discussed. 

For the hydrogen-air case, a very good agreement 
was obtained between the predicted H20 values and 
the CFD results. In addition, the predicted val- 
ues at the combustor exit plane showed very good 
agreement with the experimental value. The cost 
of doing the ILDM prediction was only a few min- 
utes to generate the look-up table and a negligible 
amount of CPU time to perform the interpolation. 
For this case then, the procedure looks very promis- 
ing. 

For the Methane case, the C02 and H20 results 
were very good, except at low concentrations. This 
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is caused by a combination of two factors. Firstly, 
the low concentration values correspond to very 
lean or unburnt compositions. Under these con- 
ditions, it is likely that mixing is a more important 
process than reaction, and so the ILDM table will 
not represent the state of the composition well. It 
is possible to include mixing in the table generation 
process, which may help this problem, but it has not 
been done in this case. The second factor is that the 
reduced mechanism may not agree with the ILDM 
scheme in certain regions. This would result in the 
compositions at the same progress variable differing 
between the two methods, and so the mapping not 
agreeing. 

The comparison of the NO values for the methane 
case does not exhibit the nice one-to-one agreement 
of the major species. However it should be con- 
sidered that the majority of data points do agree 
within an order of magnitude at the worst, and the 
high concentration values by a factor of 2. In ad- 
dition, the values obtained by the CFD results can 
by no means be considered a true solution. 

For the JetA case, there is not a very good agree- 
ment between the predicted and the CFD results for 
the major species. In this case, it is probably the 
full mechanism from which the ILDM table was ob- 
tained that is at fault. This mechanism is very old, 
and contains many outdated values for the individ- 
ual rate expressions. In addition, the values of NO 
are only obtained by a Zeldovich mechanism, rather 
than the more complete mechanisms used in the hy- 
drogen and methane cases. Despite the bad agree- 
ment though, the predicted NO values did show the 
right trend. 

Because of the non-itterative nature of this post- 
processing technique, the transport, (or advection) 
of NOx is not accounted for explicitly. To illus- 
trate this point, consider a ID combustion prob- 
lem. Initially, the fuel and air have mixed and start 
to react. Moving downstream the composition ap- 
proaches the fully burnt state, and after sufficent 
time reaches it. In this example, the NOx post- 
processor would show the correct rise in NOx level, 
until the fully burnt state was reached. After the 
fully burnt state was reached, the post processor 
would predict a constant NOx concentration. How- 
ever in reality, the NOx levels would continue to 
climb because the time scales of the NOx produc- 
tion are much slower than those of the major reac- 
tion components. The conclusion to be drawn from 
this is that this method will work best for short 
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residence time problems. 

For another case, consider the previous exam- 
ple, but at some distance downstream add a dilu- 
tion jet. For this example, the NOx post-processor 
would show the correct evolution of the NOx level 
until the dilution jet. After the dilution jet the 
postprocessor would predict a sudden drop in the 
NOx concentration, due to the sudden change in 
composition. However what one would expect to 
happen would be a gradual change in the NOx lev- 
els, due to the slow NOx reaction rates. From this 
example, it can be concluded that the post proces- 
sor can be expected to fail when there is sudden 
changes in the composition of the mixture, except 
at the initial conditions. 

Conclusions 

A very fast and reasonably accurate technique 
has been developed to post-process CFD combus- 
tion results and obtain an estimate of the NOx lev- 
els throughout the field. The method offers a signif- 
icant improvement over iterative methods in CPU 
time required, while providing a better prediction 
than equilibrium techniques. However the method 
is restricted to cases which have relatively short 
residence times, and that do not exhibit abrupt 
changes in composition or temperature. 

For the three test cases considered in this paper, 
the hydrogen-air and methane-air examples showed 
good agreement with available data. The JetA ex- 
ample did not perform as well, but the deficiencies 
can be blamed on the poor performance of the full 
mechanism used to create the look-up tables. 

Further work will be required to validate this 
method, as well as the development of better full 
mechanisms for JetA-air reaction. 
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CFD Result 



Figure 2: CFD result vs. ILDM prediction for mass fraction of C02 in CH4-air combustor. Mixture 
fraction range 0.02 to 0.07. Approximately 10,000 samples. 
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CFD Result 



Figure 3: CFD result vs. ILDM prediction for mass fraction of H20 in CH4-air combustor. Mixture 
fraction range 0.02 to 0.07. Approximately 10,000 samples. 
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ILDM Result 


Figure 4: CFD result vs. ILDM prediction for mass fraction of NO in CH4-air combustor, 
range 0.02 to 0.07. Approximately 10,000 samples. 


Mixture fraction 
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CFD Result 


.22 



Figure 5: CFD result vs. ILDM prediction for mass fraction of C02 in JetA-air combustor. Mixture 
fraction range 0.02 to 0.07. Approximately 10,000 samples. 
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CFD Result 



Figure 7: CFD result vs. ILDM prediction for mass fraction of NO in JetA-air combustor. Mixture fraction 
range 0.02 to 0.07. Approximately 10,000 samples. 
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